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We investigate the effects of an electric current on the width of a stationary reaction zone in an irre- 
versible A~ +3^^ — > C reaction-diffusion process. The ion dynamics of the electrolytes A = (^4+, ^~) 
and B = {B^,B~) is described by reaction-diffusion equations obeying local electroneutrality, and 
the stationary state is obtained by employing reservoirs of fixed electrolyte concentrations at the 
opposite ends of a finite domain. We find that the width of the reaction zone decreases when the 
current drives the reacting ions towards the reaction zone while it increases in the opposite case. The 
linear response of the width to the current is estimated by developing a phenomenological theory 
based on conservation laws and on electroneutrality. The theory is found to reproduce numerical 
solutions to a good accuracy. 
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I. INTRODUCTION 

The reaction-diffusion process A + B ^ C combined 
with the ensuing evolution of the reaction product C (e.g. 
precipitation) underlie the explanation of a large number 
of physical, chemical, and biological phenomenaii^i^iii. 
An interesting aspect of the reaction-diffusion part of 
these processes is the presence of reaction zones. They 
are formed either because the reagents A and B are ini- 
tially separated^ or because spatial inhomogeneities exist 
in the initial distribution of the reagents®. These zones 
are important since by determining where and when the 
reaction product emerges, they set the stage for the tem- 
poral and spatial evolution of Cs. Accordingly, the mo- 
tion of these fronts, the spatial distribution of the rate of 
the production of C, and the width of the reaction zones 
have been much investigated. They are known theoret- 
ically for the case of neutral reagents^i^iSiiSiii^ and the 
theories have been verified in experimentai^'i^iii^'i^ii^. 

In realistic situations, however, the reagents A and B 
are often electrolytes which dissociate, 

A-^A++A-, B^B++B-, (1) 

and the reaction takes place between the 'active' ions 
which, for definiteness, will be taken below to be A~ and 

B+: 

A-+B+^C. (2) 

Although the counterions v4+ and B^ are not reacting, 
they influence the dynamics significantly through the 
electroneutrality constraint so that the task of charac- 
terizing the front becomes much more involved^^. 

Some of the front properties such as the spatial loca- 
tion and the reaction product distribution (but not the 
width of the reaction zone) have nevertheless been ob- 
tained for ionic reactions in one-dimensional geometryi^. 
Furthermore, these studies have been extended to ions 



being driven by an electric current or by an external po- 
tential differenced^. The driven systems are of special in- 
terest when the Cs undergo phase separation, since then 
they may be used to design bulk precipitation patterns. 
Indeed, it has been shown recentl y^^i^" that a flexible 
control of precipitation patterns can be achieved through 
controlling the reaction zones by appropriately designed 
time-dependent currents. 

From the technological point of view, controlling pre- 
cipitation patterns becomes relevant if the patterns can 
be downsized to the submicron range. Since the width 
of the reaction zones is one of the limiting factors in 
downsizing, it is clear that one should understand how 
to control it. The studies of the width for neutral 
reagents^i^iii^ suggest that the parameter strongly af- 
fecting the width is the reaction rate constant. It is, how- 
ever, not a parameter we can adjust, thus other means of 
control should be found. Since electric currents turned 
out to be useful in manipulating pattern o^^i^° , it is natu- 
ral to ask if the width could also be controlled by them. 
This is the question we address in this work. 

In order to simplify the task, we restrict our study to a 
one-dimcnsional reaction-diffusion process on a finite in- 
terval. The concentrations of the electrolytes are fixed at 
the boundaries and, furthermore, a current generator is 
attached so that a constant current flows through the sys- 
tem (see FiglHand Figl2|). For this setup, we derive the 
reaction-diffusion equations in the long-time limit when 
the stationary state is reached. Solving these equations 
numerically (and, in some limits, analytically as well), we 
obtain the width of the stationary reaction front, w{J), 
as a function of the electric current, J. 

The J = case has been investigated earlier^° and we 
recover the zero current width, w(0), obtained in that 
work. For J ^ 0, our general flnding is that a forward 
current (a current that drives the reacting ions towards 
the reaction zone) reduces the width of the stationary 
reaction zone, while a current of opposite polarity (back- 
ward current) increases the width. For small currents the 
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change compared to the zero-current case is proportional 
to the current [w{J) — w{0) oc J] with a proportional- 
ity constant that can be estimated from rather general 
reasoning and the results are found to compare favor- 
ably with the numerical solutions. Although our results 
concern the width of stationary states, we present some 
arguments that the quasistationary nature of the diffu- 
sive fronts allows the derivation of the dynamics of the 
width in some time window in case of moving fronts as 
well. 

The paper is organized as follows. In Sec. II, we intro- 
duce the problem and discuss the equations describing 
the dynamics of the ions. Sec. Ill contains the deriva- 
tion of the model equations for the particular case of 
the stationary state. Analytical and approximate solu- 
tions are obtained both for symmetric setup (Sec. IV) and 
asymmetric setup (Sec. V). Possible applications of the 
stationary state results to moving reaction front are pre- 
sented in Sec. VL and conclusions are drawn in Sec. VII. 



II. THE PROBLEM AND THE MODEL 

Fig [T] displays the setup for producing a stationary re- 
action zone in the presence of an electric current. Two 
dissociating electrolytes A and B [see Eq.(IT])] are dis- 
solved in a column of gel where their transport is re- 
stricted to be diffusive. The opposite ends of the columns 
are connected to reservoirs of A and B, respectively, thus 
we fix the electrolyte concentrations oq and bo at the 
boundaries. Furthermore, we maintain a constant elec- 
tric current through the system by attaching a current 
generator to the ends of the column separated by a dis- 
tance L (in another possible setup one keeps a fixed po- 
tential difference between the ends). 




FIG. 1: (Color online) Experimental setup for generating a 
stationary reaction front in the presence of an electric current. 
The details are described in the text. 

As the reaction-diffusion dynamics proceeds, a reaction 
zone (a spatial region where the reaction A~ + 5+ C 
takes place i.e. where the rate of the production R{x, t) 
of Cs is nonzero, see Figl^]) is formed. In the long-time 
limit the system relaxes to a stationary state and the 



front becomes stationary as well, R{x,t) R{x). We 
shall be interested in the stationary state properties of 
the front, and particularly in its width w which we can 
define e.g. through the second moment of the production 
rate 



J dx R{x) 



(3) 



where the integrals are over the interval [—L/2, L/2] and 
X f is the center of the reaction zone 



J dx X R{x) 
J dx R{x) 



(4) 




FIG. 2: (Color online) Characteristic stationary profiles of 
the ion concentrations in the presence of an electric cur- 
rent. The rate of production R{x) of Cs is measured in 
units of kao where k is the rate constant of the reaction, 
and a magnification factor 10^ is used for better visibility, 
i.e. R* = 10'^-R(a::)/fcao. The position of the front is denoted 
by Xf. 

In order to calculate w, we shall use reaction-diffusion 
equations which amounts to a mean-field description of 
the problem. In addition, the following simplifying as- 
sumptions will be made: 

(i) The system is treated as one- dimensional^ i.e. all 
the relevant quantities are assumed to depend on a single 
spatial coordinate x, with —L/2 < x < L/2 being the 
spatial extent of the system. 

(ii) We consider instantaneous 100% dissociation of the 
electrolytes A and B according to Eq. ([1]) (this is the 
assumption of "ideally strong" acid and basis). Accord- 
ingly, only A^ and B^ ions (and the reaction product 
C) are present in the system. 

(iii) The dynamics of the inert reaction product C is 
assumed to have no feedback on the dynamics of the 
reagents. 

(iv) The equations are assumed to satisfy electroneu- 
trality^ meaning that the local charge density is zero on 
relevant space scales ( in particular on the scale of the 
width of the reaction zone), i.e. 



■ni{x,t) = 0, 



(5) 
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where Zi is the charge of the i-th ion - in terms of the 
elementary charge q - and ni(x,t) denotes its density. 

(v) Monovalent ions, i.e. \zi\ = 1 for all i will be con- 
sidered. 

(vi) We assume equal diffusion coefficients for the ions, 
i.e. Di = D for all i. 

(vii) The boundary conditions follow from the assump- 
tion that there are infinite reservoirs providing fixed con- 
centrations of the ions at the border . Denoting the con- 
centration profiles of the and ions by a^{x, t) and 
b^{x,t), respectively, the boundary conditions then can 
be written as 

a+{-L/2,t) = a-{-L/2,t) = aa (6) 
b+{-L/2,t) = b-{-L/2,t) ^0 (7) 
a+{L/2,t) = a-{L/2,t) ^0, (8) 
b+{L/2,t) ^b-{L/2,t) ^bo. (9) 

Since we are interested in the stationary state, the 
detailed initial conditions are irrelevant provided the 
reagents are separated initially. Some dynamical prop- 
erties, but not the width, have already been investigated 
for a similar setupi^ for the case uq <C bo where the initial 
conditions lead to a moving front. 

The general evolution equations for the concentration 
profiles a^{x,t) and b^{x,t) were derived ini^. For the 
simplified case we are considering here, they acquire the 
following form: 



da {x,t) 

m 

db+{x,t) 
di 

(?a+(x, t) 
di 

db-{x,t) 
di 



d^a- ,d(a-E) 
D-^ + X- ' 



dx^ 



dx 



ka-b+ 



DSC-x'X^-ka 



dx"^ 



dx 



d^a+ .dia+E) 

E)—- — — A- 



dx"^ 



dx 



d%- dib-E) 



dx^ 



dx 



(10) 

(11) 
(12) 

(13) 



Here D are the diffusion coefficients of the ions, fc is the 
rate constant of the reaction +-B^ ^ C, A = DF/RT 
with F = qN^ being the Faraday's constant (i.e., the 
electric charge transported by a mole of monovalent pos- 
itive ions), R is the universal gas constant, and T is the 
temperature. 

The scaled local electric field XE{x, t) in the above 
equations is obtained from the local electroneutrality as- 
sumption, and it is given byiS^J^: 



XE{x,t) 



Jit) 



q {a+ + a- +b+ + b-) 



(14) 



where J(t) is the electric current density, flowing through 
the system. In view of the electroneutrality condition, 
J{t) is divergence- free, i.e., for the one-dimensional case 
it can depend only on time. 

At this point, there are two ways to proceed. If a 
current generator is used and the current is fixed to be 



constant, J{t) = J, then Eq. (fT4|) determines E{x,t) and 
thus Eqs. ljlOliri)) together with the boundary conditions 
Eqs.lOn]) form the closed set of equation to be solved. 
This is the case we shall treat in detail below. 

Experimentally, it may be more convenient to maintain 
a constant voltage difference U = V{L) — V{0) instead 
of a constant current. In that case, integrating Eq. fHl) 
yields 



L/2 L/2 

dxE{x, t) — I dx 



J{t) 



-L/2 



-L/2 



Xq{a+ + a- +b+ + 6" 



= -U 



(15) 

and so J{t) is given through U . Then, substituting J{t) 
into Eg. p^ . one finds the scaled field XE 

/ dx {a+ +a- +b+ + b-y^ 

-L/2 

thus arriving at a closed set of Eqs. (|10m"^ . 

It should be noted that the constant U and constant J 
cases are equivalent provided only the stationary states 
are considered. Indeed, the concentrations are time in- 
dependent in the stationary state and, consequently, it 
follows from Eq.ljTS]) that constant U implies that J is 
independent of time, thus we returned to the constant 
current case. One cannot give, however, a simple rela- 
tion between U and J since, as can be seen from Eq.([T5l). 
the proportionality constant depends on the stationary 
state reached for a given U or J. 



III. STATIONARY STATE 

The equations for the stationary solution ni{x,t) = 
ni{x) are obtained by setting the time derivatives to zero 
in Eqs. (|10m3p . To make the equations more transparent 
one can introduce dimensionless quantities by measuring 
X in units of L and the concentrations in units of cq. 
Then, substituting Eq.(IT4]) with J{t) = J into Eqs. lflB - 
1131) (i.e. considering the constant current case), we obtain 



(a ) +Jo 
(6+)" - Jo 



a+ + a" + 6+ + b- 
b+ 

a+ + a- + b+ + b- 



- koa-b+ = (17) 

- koa-b+ = (18) 

u -r u J 

(a+)"-Jo( ^^ , , =0 (19) 



a+ + a- -\-b+ + b- 
b- 



^a+ + a- -\- b+ + b 

where prime denotes the spatial derivatives, and 

JL kaoL"^ 
Jo = — — , ko 



= (20) 



qDao 



D 



(21) 
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In order to solve the above equations the rescaled cur- 
rent Jq and rate constant fco, and the boundary condition 
bo/ao should be given. In principle, we could proceed 
then by numerically solving the equations. 

We would like, however, to find first some analytical 
estimates of w{Jq). For this purpose, we begin by consid- 
ering a symmetric setup where oo — bo and the ko oo 
limit is taken. In this case, Jq remains the only con- 
trol parameter, and the reaction zone becomes pointlike. 
The limit is nevertheless relevant. First, because the real 
reaction zones are narrow, in general. Second, because 
the analytical estimates of the sensitivity of concentra- 
tion profiles to small currents can be used to estimate 
the width for fcp finite but large. Once the phenomeno- 
logical estimate of w{Jq) is obtained, we carry out nu- 
merical investigations as well to judge the accuracy of 
the phenomenology. In a final step, the results will be 
generalized to the ao 7^ &o case. 



IV. SYMMETRIC SETUP (ao = 60) 

A. Point reaction zone (infinite reaction rate) 

The reaction zone is pointlike for fco ^ 00, and it fol- 
lows from the symmetry of the ao — bo setup that the 
position of the reaction zone is at Xf — 0. The con- 
centrations of the reacting ions satisfy a^{x > 0) = 
and 6+ (a; < 0) = and so we have two sets of 
equations for the two sides of the front. We begin 
by considering the left-hand side {x < 0) and obtain 
the concentration on the right-hand side by symmetry, 
namely &+>o(a;) = a-^^i-x), b-^„{x) = a+<o(-a;) and 

The equation for a; < are found by setting 6+ = in 
Eas. (117ll9l20l) : 



(a") +Jo 
(a+)" - Jo 



a 


a+ 


+ a- - 






a+ 




a+ 


+ a- - 






b- 




a+ 


+ a- - 


f 5- 



(22) 
(23) 
(24) 



The boundary conditions to the above equations at 
—L/2 [see Eqs. (|6l7l )] do not change while, at the reaction 
zone, the following boundary conditions must be used 



6-(0) 



0, 

«+(0), 



(25) 
(26) 
(27) 



The first equality follows from the infinite rate constant. 
The second one is the electroneutrality condition {a^ — 
a~ —b^ ~ 0) employed at the reaction zone. Finally, the 
third condition follows from the requirement that, for the 



counterions, not only their concentrations but also their 
derivatives should be continuous across the reaction zone. 

The first step in solving the equations for a; < is the 
application of the electroneutrality condition (a"*" — a" — 
&- = 0) to Eq.jlSl). It yields (a+)" = 0, thus resulting in 
the following solution for a+: 



a+{x) = 1 + A(2.T + 1) 



(28) 



where we used the boundary condition at a; = —L/2. 
The integration constant A remains undetermined at this 
stage. 

Having the solution for a^{x), and using the elec- 
troneutrality condition in Eq. (l22p . we find now that a~ 
satisfies the following equation 



Jo 



2 l + A{2x + l) 



= Ii 



where Ii is another integration constant. 
The general solution of equation (j29|) reads 



a- (x) = Ci [I + A{2x + 1)] + C2 [1 + A{2x + 1)]^ 
and the boundary conditions provide Ci and C2 

(l + A)^-^ 



C2 



1-{1 + A) 



A-l 



1 



1-{1 + A) 



A-l 



(29) 

(30) 

(31) 
(32) 



with A = -Jo/4A. 

Once a~ and a+ are known, b^ is obtained from the 
electroneutrality condition b~ = — . Thus what 
remained is to determine the integration constant A. We 
can find A by combining the electroneutrality condition 
with the boundary condition Eq. (|27p to arrive at 

(33) 



(«~)'L=o = 2(a+)'L ^ 

Substituting now Eqs. and ([201) into Eq 

obtain A as the solution of the following relation: 



AA' ^ 



A) 



(34) 



We have thus determined all the concentration profiles 
for X < and, as mentioned above, the profiles for a; > 
can be obtained from symmetry considerations. 



B. Reaction rate and the slope of a at a; = 

A simple question one might ask about the ko ^ 00 
case is the following: How does the rate of reaction 
changes when the current is switched on? In order to 
calculate this quantity, the slope (a^)' needs to be eval- 
uated at the reaction zone (a; = 0). As we shall see, 
the same slope will also be important in the next Sub- 
section IIV CI where a phenomenological estimate of the 
width w{J) will be carried out. 



5 



The derivative (a YI^^q can be calculated by not- 
ing from Eq.® that (a")'|^^o = 2 (a+)'|^^o and that 
(a+)' = 2 A follows from Eg. ([28|l . As a result, we arrive 
at 



(a- 



= AA. 



(35) 



The rate of reaction is the flux of the reacting ions j into 
the reaction zone. It is the number of A^ ions reaching 
X = in unit time, and is obtained as 



Joa~ 



(36) 



x=0 



where we note that j is also scaled quantity, i.e. it is 
measured in units of Dag/L. Since a"|j.^Q = 0, the 
second term on the right-hand side is zero. The first one 
is obtained from Eq. (|35p . thus yielding 



j(Jo) = ~(a-)'|_„ 



AA(Jo) 



(37) 



where A{Jo) is given by solving Eq. p4|) . 

Analytic expressions of j{Jo) can be developed for 
small Jo since we can expand the solution of Eq. ([M|) 
in powers of Jo ^ 1- To first order in Jq, we find 
yl(Jo) = -1/2+ (1 - ln2)Jo/4 and consequently 



jiJo) 



(a- 



2- (1 - In 2) Jo 



(38) 



Introducing the scaled flux jo of A~ ions in the absence 
of electric current [jo = 2), we can write the expansion 
as 



JiJo) = jo 



1 — a— 

jo 



(39) 



where a = 1 — In 2 > 0. 

Let us consider now the sign of the first order contri- 
bution. Note that Jq ~ J and it follows from Eqs. (|10lir51) 
that J < means that the negative (positive) ions are 
driven to the right (left). In this case, the reacting ions 
are driven towards the reaction zone and the first or- 
der correction is positive and, consequently, the reaction 
rate increases [j(Jo) > jo]- We shall call the J < elec- 
tric current as forward current. The reagents are driven 
away from the reaction zone in the opposite case ( J > 0, 
backward current) and, as expected, the reaction rate de- 
creases [j(Jo) < jo]- 

An important point to recognize here is that the change 
in the reaction rate does not come directly from the drift 
term in the flux of particles, it comes indirectly from the 
change in the diffusive flux i.e. from the change in the 
slope of the reaction proflles near the reaction zone. We 
shall see a similar effect when calculating the width of 
the reaction zone. 

In closing this section, we display the diffusive flux 
—D(a~)' in terms of the original variables 



-£1(0" 



\x=0 



(l-ln2 



J 



]0-OL- 



J 



(40) 



where jo = ia^D/L is the unsealed diffusive flux in the 
absence of electric current. We shall need the above ex- 
pression in the phenomenological arguments of the next 
section which would be less transparent using the scaled 
variables. 



C. Finite reaction rate: Phenomenological 
considerations and numerical solutions 

For the limiting case of an infinite rate constant (fco — > 
oo), we have analytical solutions for the concentration 
profiles for arbitrary Jq. The physically relevant case, 
however, is the one with finite fco where there is little 
hope to find exact solutions. For this case, we developed 
phenomenological considerations which have their roots 
in the success of similar arguments for simpler cases^i^, 
and we expect it to be performing well at least for small 
currents. As we shall see, the validity of the phenomeno- 
logical argument is supported by numerical integration 
of the original equations. 

The main idea is to use the balance equation for the 
reacting ions, namely to equate the number of reactions 
per unit time to the flux of reacting ions towards the 
reaction zone. To do this, one needs the values of the 
concentrations and their derivatives. They are estimated 
by assuming that the concentration profiles are smooth 
functions and, for finite but large /co, their slopes in the 
reaction zone can be approximated by the fco ~* oo val- 
ues. 

The balance equation can be obtained by integrating 
eq. (|17p through the reaction zone i.e. from —w to w 



-D- 



, da 
dx 



J 

q 



= fc'C 



b+l 



(41) 

Here the upper bars represent the spatial averages over 
the reaction zone and, in the spirit of mean-field approx- 
imation, the average of the product a^6+ has been re- 
placed by the product of averages a~ -6+. Furthermore, 
the contribution from the upper limit of the integrals 
on the left-hand side have been neglected since a~ and 
{a~y become zero at the right end of the reaction zone 
[a-~{w) « 0, (a-yiw) « 0]. 

The first term on the left hand side is the diffusive flux 
which has to be evaluated to first order in J. Our approx- 
imation consists in replacing the slope (a~)' at x = —w 
by {a~y at X = and use the expression obtained in the 
k oo limit given by eq. (|40p . The second term can be 
approximated by replacing the denominator by ao (it is 
exact in the k —^ oo and J — *■ limit). As we shall see, 
this term turns out to be negligible thus the details are 
immaterial. As a result, we find 



J J 

Jo — Oi a 

q qao 



= ka- ■ 6^ 



(42) 



The values of a and a~ ^ 6+ can be estimated by 
noting that a^{w) « thus the function at a; « is 
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approximately given as (a )'(—?«). Thus, we have 

\—w ^ ' la;— U 



^43) 

where we neglected multiplicative factors of order 2, and 
used again Eq. ((40|) for evaluating (a~)'|g. 

Substituting the above expressions into ((42|) we obtain 



Jo - a 



J 



q qaoD 



JQ 



[JO - «- 



(44) 



For zero electric current (J = 0), the solution of this 
equation is 



Wo 



kjo 



1/3 



(45) 



which is a result obtained in a study of the reaction zone 
in case of neutral reagents^°. 

Expanding now the solution of Eq.([331) to first order 
in J/qjo, one finds 



w = Wo 



1 + 




2wo^ 


-1 


(1- 


3L ) 





(46) 



Since, for the usual situation of large fc, we have wo L, 
the last term on the right-hand side (whose origin can be 
traced back to the drift term in the flux of the ions) can 
be neglected. Thus we arrive at the final form for the 
width of the reaction zone: 



w 1 - In 2 J 

— = 1 + ^^ 

Wq 3 q]o 



(47) 



This is the central result of the paper. It tells us that the 
width of the reaction zone decreases for forward currents 
( J < 0) while it increases for current of opposite polarity. 



0.8 




1 + 0.08 J/(qjg) 



X numerical result 

— linear fit 

— theoretical approximation 



J/(qjo) 



1 



FIG. 3: (Color online) The small-current expansion of the 
width from numerical integration of the Eas. (l22l24p is com- 
pared to the phenomenological theory of Subsection IIV CI 
The width is scaled by its value at J = while the current is 
scaled by gjo where jo is the particle flux at J = 0. 

The numerical factor in front of J/qjo is of course a 
question and this is why we carried out the numerical 



integration of the differential equations. The equations 
are simple and the numerics does not pose any problem. 
FigE] shows the comparison of the numerical results with 
the phenomenological theory. The agreement is surpris- 
ingly good in view of the simplicity of the phenomenolog- 
ical arguments we used. Thus our assumptions entering 
the derivation of w{J) in the small J and large k limits 
appear to be justified. 



V. ASYMMETRIC SETUP 

We would like to control the width and so, it is useful 
to consider a more general case provided by asymmetric 
boundary condition oq 7^ bo since it introduces a new 
control parameter 



Q = 



ao 



(48) 



The study of this case follows the steps of the symmetric 
case. First, the point like reaction zone (fco — > 00) is 
solved exactly, then the phenomenological consideration 
are repeated, and finally, the numerical simulations are 
used to check the validity of the phenomenology. 

As before, for ko ^ 00, the left and the right sides of 
the front can be treated separately. For the left hand side 
{x < Xf) one can use the equations (|^ . (|^ and pi)) 
while for the right hand site {x > Xf) the equations are 
as follows 



b+ 



(a+)" - Jo 
(b-)" +Jo 





f 6- 


+ a+ 




a+ 






hb- 


+ a+ 




b- 




U+4 


-b- - 





= 
= 
= 



(49) 
(50) 
(51) 



and the electroneutrality condition for x > Xf reads 

a+ + b+ -b- =0 . (52) 

These equations are solved with the following bound- 
ary conditions for the dimensionless quantities. At the 
outer boundaries we have a^(— 1/2) = 1, &~(— 1/2) = 0, 
a+(l/2) = 0, and &^(l/2) = Q. At the front {x = xj), 
one has a~{xf) = 0, b'^{xf) = and a^{xf) — b^{xf). 
Finally, the conditions for the slopes at x — Xf are: 



(53) 
(54) 
(55) 

"J- " t 

where xj (x^ ) denotes the left (right) side of the reaction 



{a-y 






(«+)' 




= 


{b-y 


Xj 





zone. 



Since the solution of the above problem is similar to 
the symmetric case (though the algebra is much more 
tedious), we present only the main results. 
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Due to asymmetry, the position of the front Xf shifts 
from X ~ and one finds that it is given by 



l-Q 
2(1 + 0) 

0-0'- 



Jo 



20(1 + 0)' 

1 



(56) 



In 



1 + 



0' In 







1 + 



An important quantity for the phenomenological argu- 
ments is the slope of the concentration of the reacting 
ion on the left-hand side of x/ . It is found to be 



{a-)\^^ = -{l + Q) + a{Q)Jo 



with 

a(0) = 1 



■In 



20 \l + Q 







In 







1 



(57) 



(58) 



Fig. m shows the dependence of a on and one can ob- 
serve that a{Q) has a minimum at = 1 in the 
symmetric setup. As expected the expression (|58p is in- 
variant under the exchange of the boundary concentra- 
tion values, ao and bo, i.e. a{Q) = a{l/Q). 



0.32 



0.315 



0.31 



l-ln2 



0.305 



FIG. 4: (Color online) Q dependence of the parameter a. 

For finite reaction rate, one can follow the steps of 
the phenomenological arguments of the symmetric case, 
and calculate the first order expansion oiwl J) in J/{qjo) 
where jo is given now by 




jo = (ao + k))D/L . 



The final result is 



^ ^ 1 + "(Q) J 

Wo 3 qjo 



(59) 



(60) 



suggesting that the width of the reaction zone is more 
effectively reduced when the system is asymmetric. 

The predictions of the phenomenological arguments 
can again be compared with the results of a numerical 
integration of the differential equations. The agreement 
is qualitatively similar to the symmetric case and the ten- 
dency of an increasing impact of the imposed current on 
the reaction front width when enhancing the asymmetry 
has been verified. 



VI. CONSIDERATIONS ABOUT THE WIDTH 
OF A MOVING FRONT 

Although our motivation for this work comes from con- 
trolling precipitation patterns which emerge in the wake 
of a moving reaction-diffusion fron-feiSiSS, as a first step, 
we considered the controllability of the width in a sim- 
pler stationary state. Clearly, it should be now clarified 
to what extent the results are generalizable to the case 
of moving fronts. 

The first step in the generalization is the realization 
that diffusion fronts slow down with time (their velocity 
is proportional to l/\/t). This means that the system 
enters a quasistationary regime where stationary state 
considerations usually yield correct result. As a rele- 
vant example here, let us examine the time evolution 
of the width of the reaction zone wo{t) emerging when 
the reagents are initially separated and no current is 
present [a{x < 0,i = 0) = ao, h{x < 0, i = 0) = 0, 
a{x > 0, i 0) = 0, 6(x > 0, i = 0) = 6o < ao, L — > oo]. 
In this case, the front moves diffusively and the width 
grows with time as- WQ{t) ~ t^/^. This time-evolution 
can be derived easily from the stationary state resulti^ 
w ~ {D"^ /kjoy/^ [see ea. ps]) ] by the following argument. 
The width of depletion zone (the region where the con- 
centration profiles are significantly different from their 
initial values) grows with time as Wd{t) ~ Vot. Then 
the diffusive current towards the reaction zone is obtained 
as jo(t) ~ D{ao + foo)/Wrf(i) and one finds the 

desired time-evolution as 



Wait) 



kjo{t) 



1/3 



,1/6 



(61) 



The second property of the front needed for the gener- 
alization is that the motion of the front remains diffusive 
even if the ionic nature of the reagents are taken into 
account and a small current is switched on^^. Assuming 
that the diffusive nature of the front implies quasista- 
tionarity, one can go through the estimates of the various 
terms in the reaction-diffusion equations provided those 
terms are restricted to the neighborhood (the region of 
depletion zone) of the reaction front. Using the quasista- 
tionary time evolution of the diffusive current Jq ^ 
one finds then that the neglected time derivatives are in- 
deed small provided w{t) <C Wdit). Furthermore, follow- 
ing along the same steps which lead to the estimates of 
w/wq in the small electric current limit [see the deriva- 
tions of eqs.lgTl) and ((101) in Sec UVCl and Sec|V], one 
finds 



w{t) « Wo{t) 



1 



J 



Qjoit) 



(62) 



where a is a numerical constant of the order one. An 
important conclusion to be drawn from the above ex- 
pression is that the direction of the current has the same 
effect as before. Namely, the forward current decreases 
the width while the backward current increases it. 
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Unfortunately, the decrease of jo ~ t^^^^ makes it clear 
that the expansion of w in J/qjo cannot be valid for large 
times. It is valid only at intermediate times when the 
quasistationary regime has already set in but the electric 
current is still smaller than the diffusive current. 

Clearly, we cannot claim here that the width of a mov- 
ing front in the presence of a current has been under- 
stood. We believe, however, that the solution should 
have the properties described above in a well-defined 
intermediate-time regime. Finally, we should also note 
that, in the control of precipitation patterns^, one em- 
ploys actually time-dependent currents and one has some 
freedom in the design of those current. Consequently, 
the effects of variable currents on the width, and the 
use of the freedom in the choice of current to control 
the width provide additional challenges which should be 
treated once we arrived at basic understanding of the 
simpler problems. 

VII. DISCUSSION AND FINAL REMARKS 

We have studied the effects of an external electric cur- 
rent (or potential difference) on the width of the reaction 
zone in an irreversible A~ + — s- C reaction-diffusion 
process. The aim was to find a flexible control of the 
width since this would help in the attempts of down- 
scaling the patterns generated by reaction-diffusion pro- 
cesses. 

Unfortunately, we have to conclude that the width of 
the reaction zone is not very sensitive to the presence of 
an externally imposed electric current or potential differ- 
ence. In the symmetric setup for a typical values J/qjo 
of order —1, the reduction of the width compared to the 
zero-current case is of the order of 10% at best (see Figl3]) 
and, furthermore, even a strong asymmetry in the bound- 
ary conditions (Q — > oo) allows only a reduction of the 
order of 20% [see Eqs. ipiSDl) ]. 

Clearly, the new attempts at control should start now 
with the reexamination of wq. In order to see the prob- 
lem, let us write Ea. (|l5)) in terms of the characteristic 
concentration and lengthscale of the system 

^,^(£lY' (63) 



The rate constant k and the diffusion constant D can- 
not be controlled effectively. So, what remains to con- 
sider is the gradient of concentrations, ao/L, and one 
can see from Eq. (p5)) . that small width is obtained by 
large concentration gradients. Actually, the wet stamp- 
ing metho d^^'^^ employs large gradients (ao ~ lA/ and 
L « 100/i) to achieve patterns with features at submi- 
cron scales. The problem with wet stamping, however, is 
that the patterns are generated by skillfully preparing the 
boundary conditions and thus the method runs into dif- 
ficulties with the task of building three-dimensional pat- 
terns. Nevertheless, the wet stamping points towards a 
possible solution of our problem. Namely, we would have 
to design reaction-diffusion processes which can support 
large concentration gradients in a reaction zone whose 
motion can be controlled. Of course, this appears to be 
a nontrivial task. 

Although the diffusion constant D cannot be con- 
trolled, one can still think about considering materials 
where D is significantly different from those of the ions 
in gels {Dgei « 10~*m^/s). In solid state systems, like for 
example glasses, the diffusion coefficients differ strongly 
from the typical values in gelatine materials (typically 
Dsoiids ~ 10~^^m^ / s)^. Accordingly, the zero-current 
width can be extremely small, since w{0) ~ D^/-^ (see 
Eq. (61)). 

The study of glassy systems is also interesting since 
the dynamics of reactants may become subdiffusive in a 
glassy matrix. For subdiffusive systems we have D — )■ 
and, as recent studies have shown, the subdiffusive 
reaction-diffusion dynamics leads to striking changes in 
the concentration profiles of the reacting ions^i. Thus, 
combining subdiffusivity with the use of electric currents 
appears to be a promising way of producing controlled 
reaction zones of very narrow (submicron) width. 
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